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Abstract 

This paper describes the development and validation of 
an optimization procedure to systematically place tun- 
ing masses along a rotor blade span to minimize vibratory 
loads. The masses and their corresponding locations 
are the design variables that are manipulated to reduce 
harmonics of hub shear for a four-bladed rotor system 
without adding a large mass penalty. The procedure 
incorporates a comprehensive helicopter analysis to 
calculate the airloads. Predicting changes in airloads 
due to changes in design variables is an important 
feature of this research. The procedure was applied to 
a one-sixth, Mach-scaled rotor blade model to place 
three masses and then again to place 
six masses. In both cases the added mass was able to 
achieve significant reductions in the hub shear. In 
addition, the procedure was applied to place a single 
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mass of fixed value on a blade model to reduce the hub 
shear for three flight conditions. The analytical results 
were compared to experimental data from a wind tunnel 
test performed in the Langley Transonic Dynamics 
Tunnel. The correlation of the mass location was good 
and the trend of the mass location with respect to flight 
speed was predicted fairly well. However, it was noted 
that the analysis was not entirely successful at predicting 
the absolute magnitudes of the fixed-system loads. 

Nomenclature 

C T coefficient of thrust 

f objective function 

I number of constrained frequencies 

J number of tuning masses 

K number of harmonics of shear included in the 
objective function 
Mj jth tuning mass 

NP frequency or loading at N times the rotational 
speed of blade 

Sk amplitude of kth harmonic of shear 

x j location of jth tuning mass 

P k scalar design variables appearing in the ob- 



jective function and constraints 
\i advance ratio (ratio of forward flight speed to 
tip speed) 

ith natural frequency 
(Oj COj / ft 

co u lower bound on CDj 

(O ui upper bound on COj 

q rotational speed of rotor blade 

Introduction 

Since helicopter vibration is transmitted from the 
rotor blade to the fuselage through a time dependent 
shear force at the hub, methods for reducing vibration 
through reduction of hub shear have long been a 
subject of study. An example of this is vibration reduc- 
tion of rotor blades through passive control. For in- 
stance, pendulum absorbers 1 , active isolation devices 

o 2_<s 

, additional damping , and vibration absorbers which 
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create anti-resonances have all shown promise in 
reducing blade vibratory response. Historically, fre- 
quency placement has been the principal technique for 

£ 

passively reducing blade vibration . Another form of 
passive control is to alter the mass and stiffness distri- 
butions of the blade. These modifications tailor the 
mode shapes to achieve orthogonality to the airfoading 
thereby reducing the generalized force and response of 
the blade 7 . This is generally done in a late stage of the 
design process. 

The current trend in engineering design of aircraft is 
to incorporate critical requirements from all pertinent 
disciplines into an early phase of the design process to 
avoid costly modifications after a problem has been 
detected 8 . In the preliminary stage of design, a large 
number of design variables are free to be chosen in 
order to satisfy important multidisciplinary consider- 
ations. It is in this stage that passive control of vibration 
can play an important role. The design process is a 


labor-intensive effort, however, mathematical optimiza- 
tion techniques allowfor efficient and thorough searches 
of the design possibilities while satisfying a large num- 
ber of conflicting design requirements from many differ- 
ent disciplines. For example, reference 9 used optimi- 
zation methods to study the interaction of structural 
properties with airload distributions in designing blades 
for low vibration. The structural properties included 
mass and stiffness distributions. The airload distribu- 
tions included higher harmonic lift components and 
aerodynamic pitching moments which are the primary 
sources of vibration in helicopter rotor blades. Com- 
parison of the vibration characteristics from three ana- 
lytical design strategies showed the benefits of using an 
automated structural optimization procedure with a 
coupled aeroelastic analysis. Another example where 
optimization was successfully applied to the design of 
low vibration rotor blades was reported in reference 1 0 
and 1 1 where several alternative optimization formula- 
tions were investigated and their benefits revealed. 
References'! 0 and 1 1 did not use computed airloads in 
the analyses. Reference 12 discusses an optimization 
procedure for designing a low vibration rotor blade. 
Wind tunnel tests of the blade showed that the design 
proved to be better than a rotor designed using the 
traditional approach of frequency placement. A com- 
parison between the analytical results and test data 
revealed that the trends and reductions in load levels 
were predicted well but the absolute values of the loads 
at given airspeeds were predicted less accurately. 

Reference 13 described a procedure for placing 
and sizing tuning masses at strategic locations along 
the blade span to tailor the mode shapes. This proce- 
dure used formal mathematical programming tech- 
niques in conjunction with a finite element program to 
model a simplified blade and calculate the dynamic 
response. The airloads used in the analysis repre- 
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sented a set of harmonics typical of a four-bladed rotor 
system 7 . The loads did not vary with changes in the 
masses or their locations. The purpose of this paper is 
to describe the enhancement and validation of the 
method described in reference 1 3. The enhancements 
include the incorporation of a comprehensive helicopter 

14 

analysis CAMRAD/JA into the optimization proce- 
dure which yields a more realistic blade model and 
calculated airloads. The validation is accomplished by 
comparing the analytical results with experimental data 
from wind tunnel tests. 

Problem Definition and Formulation 



Figure 1 . Design variable definition for optimizing 
magnitudes and locations of tuning masses 


ables, 3k which also appear in the constraints as 
“upper limits” on the shear harmonic amplitudes, S k . 

S k /p k -1^0 k = (2) 

K represents the number of shear harmonics to be 
included in the objective function. By convention a 
constraint is satisfied if its value is less than or equal to 
zero. Consequently, the optimizer will tend to decrease 
the values of 3k to minimize the objective function but 
will also tend to increase the values of 3k to satisfy the 
constraints. This results in a compromise on the values 
of 3k which forces a reduction in the values of Sk thus 

reducing the hub shear harmonics while incurring the 
smallest possible mass penalty- Additional constraints 
include upper and lower bounds on the natural frequen- 
cies of the blade to avoid resonance as shown in 
equation (3). 


1-tOj / co H <0 


The design goal is to find the optimum combina- 
tion of masses and their locations (Fig. 1) to reduce the 
vertical hub shear. The method entails formulating and 
solving an optimization problem in which the tuning 
masses, M’s and corresponding locations, X’s are de- 
sign variables that are manipulated to minimize the 
objective function. Equation (1) defines the objective 
function, f which is a combination of vertical hub shear 
and added mass. 


f = 


fi+ ipk 


V k=l 






(D 


The objective function includes additional design vari- 


where C0j = (Oj I Cl, C0j is the ith natural frequency, 
co ui and cOjj are the upper and lower bounds on 0); 
respectively and I is the number of constrained frequen- 
cies. 

Analyses 

The analyses that are used in the procedure are the 
comprehensive helicopter analysis code, CAMRAD/ 
JA 14 , the optimization code, CONMIN 1 ® and an ap- 
proximate analysis to reduce the number of CAMRAD/ 
JA analyses during the iterative process. CAMRAD/JA 
calculates rotor performance, loads, vibration mode 
shapes and frequencies, aeroelastic stability and re- 
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sponse. In this study, CAMRAD/JA was used to calcu- 
late frequencies, airloads and hub loads. The structural 
model of the rotor is based on engineering beam theory 
for rotating wings with large pitch and pretwist. The 
frequencies and mode shapes are computed using a 
modified Galerkin analysis. The rotor aerodynamic 
model is based on lifting-line theory and uses steady 
two-dimensional airfoil characteristics provided in tables 
of section lift, drag, and pitching moment versus Mach 
number and angle of attack. The analysis includes 
unsteady aerodynamic forces from thin airfoil theory 
and the induced velocity is obtained from either a 
uniform inflow model or a vortex wake model. A detailed 
description of the theory is given in reference 14. 

CONMIN is a general purpose optimization pro- 
gram that uses the method of feasible directions for 
constrained function minimization and the conjugate 
direction method of Fletcher and Reeves for 
unconstrained minimization problems. The approxi- 
mate analysis uses a linear Taylor Series expansion to 
approximate the objective function and constraints for 
the iterative portion of the optimization procedure to 
save computational effort. 

Organization of the Procedure 

A flowchart of the optimization procedure is illus- 
trated in figure 2. The overall procedure consists of two 
nested loops. Each pass through the outer loop is 
referred to as a cycle which involves a full analysis and 
a sensitivity calculation. The sensitivity analysis in- 
cludes calculating finite-difference derivatives of the 
objective function and the constraints with respect to 
the design variables. The first step is to generate the 
model of the blade, excluding tuning masses. The 
design variables (masses and locations) then deter- 
mine where and how much mass should be placed on 
the blade. Next the modal analysis is performed and the 



Figure 2. Flow chart for optimization procedure 


airloads and hub shears are computed using CAM RAD/ 
JA, and the objective function and constraints are 
calculated. The inner loop consists of the optimization 
program CONMIN and the approximate analysis for 
calculating values of the objective function and con- 
straints. Once the inner loop has converged, the next 
cycle begins, using updated values of the design vari- 
ables. The process continues until convergence of the 
outer loop is achieved. The major improvement in the 
procedure overthat of reference 13 occurs in the outer 
loop where the airloads are calculated by CAM R AD/J A. 
The use of CAM R AD/J A enables the change in airloads 
due to changes in the design variables to be taken into 
account. 

Demonstrati on of the Method 

The model used to demonstrate the procedure is a 
4-bladed, one-sixth, Mach-scaled representation of a 
design intended to satisfy the requirements for the 
“growth” version of the U.S. Army’s UH-60A (Black 
Hawk) helicopter. Each blade (shown in figure 3) 
weighs about 3 lbs and has three sets of advanced 


4 






RCtfyiQ i- RC(3}~08 
•Mol / airfoil 


100 M7< J 5 J 5 BC(4).io»iitei > 0 -<- 

"tH i 

rr 

i . v -i—-*- 

\T\ 

> 

1.798 

VnidvMT i 
p<**T 



^ 51*413 “ 1 1 ^1 

< 5*224 


Canter 

or 

rotation 


Figure 3. Rotor blade geometry 

(Dimensions are in inches) 

airfoils, RC(4)-1 0 1 6 , RC(3)-1 0 1 7 , and RC(3)-08 1 7 . The 
planform is tapered with a -16 degree lineartwist. The 
calculations were performed in CAMRAD/JA using five 
flap/lag bending modes and two torsion modes. Each 
blade was modeled with 1 8 aerodynamic segments and 
50 structural segments. The chordwise center of grav- 
ity, aerodynamic center, and elastic axis were coinci- 
dent and located at the quarter-chord. The rotor was 
trimmed to prescribed values of thrust and zero flapping 
angles using nonuniform inflow with a prescribed wake 
geometry. The 4P blade vertical shear* is generally the 
primary source of vertical loads in a four-bladed rotor 
configuration and the 3P and 5P blade vertical shears 
contribute to the hub moments. Therefore, these three 
quantities were the major focus for reduction in this 
optimization study. The flight condition was forward 
flight at an advance ratio of jl = 0.35 and a thrust 
condition of Cy = 0.0081 which signifies a full scale 
gross weight of 18,500 lbs. 

Table 1 shows the results obtained from applying 
the optimization procedure to place three tuning masses 
along the span of the blade model. Starting with the 
baseline blade (no added mass), the procedure was 
able to reduce the 3rd, 4th, and 5th harmonics of shear 
by 8%, 8%, and 4% respectively by adding 0.338 Ibm of 

* II is customary in rotorcruft dynamics to use the 
notation N/rcv or NP to denote frequencies or 

harmonics of loading at N times the rotational speed 
of the blade. 


Table 1 . Comparison of baseline and 
optimized designs from 3 -mass and 6-mass 
optimization procedures 



3-MASB H 

6-MASS 

MAm <>!»*> £ 

0.00 

0.028 

0.059 


0.00 

0.155 

0.022 


0.00 

0.155 

0.302 

: 

0.00 

0.00 

0.304 

■ MAMXVhmy' 

0.00 

0.00 

0.280 


0.00 

0.00 

0.256 

'• mjcatwn 

Va&my 

— 

26.42 

26.80 

t^OCATIQN 

- 

23.61 

26.32 

. ixxzAtton* 

- 

29.23 

25.42 


— 

— 

27.14 


— 


26.29 

LOCATION * 

- 

— 

26.90 


18.30 

16.80 

13.90 


11.50 

ia60 

7.57 


4.80 

4.60 

3.25 


tuning mass which is about 1 1 .5% of the nominal blade 
weight. The mass was added between 42 and 52 
percent of the blade span. The baseline blade was 
originally designed for low vibration so these reductions 
in shear from the baseline design although seemingly 
modest, are considered to be significant. 

As a matter of interest the procedure was also 
applied to place six masses and results are also shown 
in table 1 . In this case, the optimizer placed all masses 
between 45 and 48 percent of the blade span. The 
procedure reduced the 3rd, 4th, and 5th harmonics of 
shear from the baseline values by 24%, 34%, and 32% 
respectively with a total addition of 1.2 Ibm of tuning 
mass. This is a sizable reduction in shear, (approxi- 
mately four times the reduction in the 3-mass case) 
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however, the added mass represents about 30 percent 
of the mass of the baseline blade. Consequently, this 
solution to the vibration problem would probably not be 
adopted in a practical design situation. Nonetheless, 
this result verifies our intuition that use of additional 
mass design variables improves vibration reduction but 
with a significant increase in the weight penalty. 


Comparison with Test Data 

Reference 1 8 describes results of wind tunnel tests 
performed in the NASA Langley Transonic Dynamics 
Tunnel (TDT) on the blade test article shown infigure 4. 
The test article, based on the model described in the 



Figure 4. Model rotor blade components 

previou s section , was designed to study passive means 
for minimizing fixed-system loads and has the capability 
of adding nonstructural masses at various intervals of 
the blade span. This test provided an opportunity to 
validate the present optimization procedure. A descrip- 
tion of the model, test set-up, and the reduced data are 
given in reference 1 8. The test involved placing a single 
mass of fixed value (.27 Ibm) at various locations along 
the blade span to determine the effectiveness in reduc- 
ing 4P hub shears and moments for several different 


flight conditions. 

Reference 18 shows several data plots of 4P nor- 
mal force as a function of mass location on the blade. 
Many of the curves are very flat and suggest that the 
mass locationdoes not significantly affectthe hub shear 
for all flight conditions. Since optimization works best 
for problems with well defined minimums, only selected 
cases were included in this study. 

The optimization procedure was applied to this 
model for the placement of a single mass at advance 
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Figure 5. Comparison of optimum location of single 
mass with test data 

ratios of 0.25, 0.30, and 0.35 and a thrust condition of 
Cj = 0.0081 . Comparisons between the optimization 
results and the test data are summarized in Figure 5. 
The measured values of the optimum locations are 
shown as 1 0 percent ranges of the blade span since the 
data was only available at 1 0 percent increments of the 
blade span. For the 0.35 advance ratio case, the 
optimization procedure predicted an optimum location 
within the range of the test data. The other two cases 
(0.25,0.30) were 11% and 12% respectively below the 
lower bound of the range. This is fairly good agreement 
considering the well-known difficulty of predicting fixed- 
system loads. 

In order to verify the mechanism by which the hub 
loads are decreased figures 6 and 7 illustrate examples 


J Test data 
-0- Analysis 
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of the changes in the calculated mode shapes and 
airloads respectively from the baseline to the optimized 
design. The simultaneous changes in the mode shapes 
and airloads result in a reduction of the generalized 
force and subsequently the hub shear. The sizeable 


Second 

bending 

mode 

shape 



Figure 6. Comparison of mode shapes from 
baseline and optimum designs 



Figure 7. Comparison of 4P airloads from 
baseline and optimum designs 


change in the airload distribution resulting from the 
changes in the design variables suggests that neglect- 
ing this effect would be erroneous. 


Table 2. Predicted versus measured results for single 
mass placement to reduce hub shear 



Advance ntio=0.25 

Advance ndo=O30 

Advarce 

ratio=Q.35 


Predicted 

Measured* 

Predicted: 

Measured*: 

mm 

Measured* 

Optimum 
location of 
mm, i/R 

0.40 

0.48 

0.44 

0.52 

056 

052 

S4 ratio 
final/mitii! 

0.W 

0.93 

0.92 

0.90 

0.90 

0.92 

$4 (IbQ 
Initial design 
mus^Q lb® : 

1020 

20.10 

7.10 

37.50 

11.15 

49.00 

S40bf) 
final design 
ms&)27 
Ibm 

9.60 

18 JO 

657 

33.60 

10.02 

45.00 


* Based on faired curve through experimental data poinu 


optimum location was in fair agreement with the test 
data and the relative change in hub shear from the 
baseline to final design was also comparable between 
the optimization and test measurements for all three 
cases. However, the absolute values of the hub shear 
were grossly underpredicted. Unfortunately, this repre- 
sents the current state of the art in prediction of fixed 
system loads in rotorcraft systems (see for example 
Refs 19,20). Although it is outside the scope of this 
paper to either investigate the reasons for or to recom- 
mend corrective action forthis shortcoming, the authors 
feel obliged to point out its existence. Further we are 
aware of a number of efforts in progress to improve the 
fidelity of fixed system loads prediction. For example, 
some of the more fruitful approaches involve the inves- 
tigations of accurate calculation of rotor airloads (Refs 
21,19,22). 


Concluding Remarks 


Additional De tails of Design Comparisons 

Table 2 shows some additional details of the three 
results for the single mass placement optimization. The 


This paper described the development and valida- 
tion of an optimization procedure to systematically 
place tuning masses along a rotor blade span to mini- 


7 



mize vibratory response. The masses and their corre- 
sponding locations were the design variables that were 
manipulated to reduce harmonics of hub shear for a 
four-bladed rotor system without adding a large mass 
penalty. The procedure incorporated a comprehensive 
helicopter analysis to calculate the airloads. This pro- 
cedure enabled the changes in airloads due to the 
changes in design variables to be taken into account. 

The procedure was first applied to a one-sixth, 
Mach-scaled rotor blade model for two cases: (1) the 
placement of three masses and (2) the placement of six 
masses. The optimized three-mass configuration re- 
duced the 3rd, 4th, and 5th harmonics of shearbetween 
4 and 8 percent. The optimized six-mass configuration 
reduced these shear harmonics between 24 and 34 
percent, although four times the mass was added over 
a smaller region of the blade. 

The optimizer was then compared to test data for 
placing a single mass of fixed value on a blade model to 
reduce the 4P hub shear for three flight conditions. A 
wind tunnel test was performed in the Langley Tran- 
sonic Dynamics Tunnel (TDT) on a blade test article. 
The analytical results were compared to the experimen- 
tal data and the trend of the mass location with respect 
to flight speed was predicted fairly well. At the same 
time it was noted that the analysis used was not entirely 
successful at predicting the absolute magnitudes of 
fixed-system loads. 
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three masses and then again to place six masses. In both cases the added mass was 
able to achieve significant reductions in the hub shear. 

In addition, the procedure was applied to place a single mass of fixed value on a 
blade model to reduce the 4P hub shear for three flight conditions. The analytical 
results were compared to experimental data from wind tunnel tests performed in the 
Langley Transonic Dynamics Tunnel (TDT) . The correlation of the mass location was 
good and the trend of the mass location with respect to flight speed was predicted 
fairly well. ^ _____ _ ___ 
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